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RESUME 


On  a  concu  un  modfele  math^matique  de  reticule  de  t&te  chercheuse 
a  modulation  d'amplitude  et  d^termin#  sa  position  d'6ciuilibre  stable  en 
presence  d'une  cible  S  double  points,  ou  &  simple  point  accompagnde 
d'une  cible  dtendue.  Le  comportement  de  ce  reticule  constitue  l'objet 
d'une  simulation  sur  ordinateur  num^rique.  (NC) 


ABSTRACT 


A  generic  AM-reticle  seeker  was  modelled,  and  its  steady-state 
equilibrium  position,  in  the  presence  of  either  two  point  tarqets  or  a 
point  and  an  extended  target,  determined  by  digital  computer 
simulation.  (U) 
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1.0  INTRODUCTION 


An  infrared  seeking  missile  employs  the  radiation  from  the  target 
as  its  homing  source;  its  gimballed  seeker  head  nrecesses  to  follow  this 
source.  Although  there  has  been  a  rapid  development  of  charge-couoled 
devices  suitable  for  processing  target  radiation,  current  threats  are 
more  likely  to  use  mechanical  scanning.  The  mechanical  scanners,  which 
have  fields  of  view  of  a  little  more  than  a  degree  to  several  degrees, 
are  either  of  the  rosette  or  reticle  types.  The  rosette-type  scans  with 
a  very  narrow  instantaneous  field  of  view,  thus  permitting  the 
recognition  of  individual  characteristics.  The  reticle  scanners  include 
those  which  provide  pulse-position  modulation,  freouency  modulation  or 
amplitude  modulation  to  encode  the  angular  position  of  the  target 
relative  to  the  seeker  axis.  In  this  document  a  seeker  employing  an 
amplitude-modulating  (AM)  reticle  is  examined;  related  reticles  are 
employed  in  certain  American  and  Soviet  missiles. 

With  two  targets  within  its  field  of  view,  the  seeker  of  a 
missile  employing  an  AM  reticle  has  a  tendency  to  point  nearer  to  the 
more  intense  one,  a  property  which  is  employed  to  decrease  the 
effectiveness  of  the  infrared  seeking  missile;  a  taraet  can  dispense  a 
decoy  flare,  which,  if  it  is  sufficiently  intense,  may  be  accented  bv 
the  missile  seeker  as  a  preferential  source  of  radiation  to  the  target 
itself.  Analyses  of  infrared  decoy  flare  effectiveness  reguire  the 
determination  of  the  position  to  which  the  seeker  of  the  attackino 
missile  will  point  relative  to  the  tarnet  and  flare  positions.  Tn 
simplify  the  calculations,  it  is  often  assumed  that  the  seeker  will 
point  to  the  weighted  mean,  that  is:  the  product  of  imaoe  intensity  and 
distance  from  the  centre  of  the  reticle  will  in  the  steady  state  be 
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equal  for  the  target  and  the  flare.  This  assumption  may  be  valid  when 
two  infrared  sources  are  far  from  the  missile,  and  the  imaqes  on  the 
reticle  are  thus  small.  But  it  will  no  longer  be  true  as  the  missile 
approaches  the  targets,  and  the  imaoes  occupy  a  substantial  portion  of 
the  reticle. 

A  model  of  a  generic  AH  seeker  is  developed,  and  its  performance 
examined  for  certain  specific  cases  which  are  tractable  by  exact 
mathematical  solutions.  As  well  as  introducing  the  reader  to  the  seeker 
behaviour,  these  examples  serve  to  validate  the  computer  programme.  A 
digital  programme  is  developed  from  the  model,  and  employed  to  determine 
the  seeker  eauilibrium  in  the  presence  of  two  point  taroets  as  a 
function  of  target  separation  and  relative  target  intensities.  This  is 
then  extended  to  include  the  case  where  one  target  is  an  extended 
source.  The  results  nay  be  of  assistance  in  the  modelling  reouired  for 
flare-effectiveness  studies. 

This  report  covers  the  work  performed  at  DREV  between  1  October 
1979  and  1  August  1980  under  PCN  21013,  "Multiple  Target  Tracking  with  a 
Reticle  Seeker”. 
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2.0  MODEL 


2.1  General 


A  generalized  schematic  of  an  infrared  seeker  is  shown  in  Fin.  I. 
The  infrared  emission  from  the  target  is  focussed  onto  a  reticle 
sninning  about  its  centre  at  a  freciuency  of  f  .  The  reticle  consists  of 
alternate  opaque  and  transparent  sectors,  which  chop  the  incident 
radiation  to  produce  a  carrier  of  frequency  p.fg  a*-  the  photocell 
output,  where  p  is  the  number  of  pairs  of  opaaue  and  transparent 
sectors.  The  chopping  efficiency  is  maximum  for  an  image  capable  of 
being  instantaneously  masked  or  unmasked  totally;  such  image  can  be 
either  a  point  or  a  straight  line  (or  its  extension)  passion  throunh  the 
centre  of  the  reticle.  This  waveform  is  modulated  at  the  spin  frequency 
by  a  semi-transparent  phasing  sector.  An  automatic  qain  control  (AGC) 
keeps  the  amplitude  of  the  carrier  constant  and  independent  of  the 
target  intensity,  which  is  a  requirement  for  linear  control.  This  signal 
is  demodulated  and  produces  quadrature  components  at  the  spin  freauencv, 
the  amplitudes  of  which  are  a  measure  of  the  angular  error  between  the 
look  angle  of  the  seeker  and  a  line  to  the  target.  These  error  signals 
are  applied,  via  torouinq  motors,  to  nrecess  the  reticle  assemblv  to 
reduce  the  error  angle  to  a  small  value. 

2 . 2  Reticle 

The  reticle  used  in  this  model  is  shown  in  Fin.  2.  It  possesses 
twelve  pairs  of  transparent  and  opaque  spokes,  and  the  inner  boundaries 
of  the  phasing  sector  are  delimited  by  semi-circles,  centred  on  opnosite 


sides  of  a  diameter  of  the  reticle,  the  radiuses  of  which  are  half  that 
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PHOTOCELL 
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FIGURE  1  -  Control  system 
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of  the  reticle.  The  phasinq  sector  has  a  transmission  of  50‘V,  so  that 
the  averaqe  siqnal  will  be  the  same  durinq  imaqe  dwell  nn  the  nhasinn 
sector  and  imaqe  dwell  on  the  spokes;  this  will  tend  to  improve 
discrimination  aqainst  large  imaqes.  Harmonic  reference  siqnals  are 
qenerated  at  the  spin  freouencv  of  the  reticle,  which  rotates  in  a 
clockwise  direction  at  a  rate  of  u>o  rad/s. 

To  illustrate  the  manner  in  which  the  reticle  encodino  is  used  to 
produce  error  siqnals,  an  example  is  now  oiven  in  which,  for  ease  of 
comprehension,  the  signal  processing  is  greatly  simplified.  Let  the 
target  image  be  at  a  distance  r  from  the  centre  of  the  reticle.  Then, 
from  the  reticle  geometry  qiven  in  Fig.  5,  it  can  be  seen  that  the  imaqe 
will  be  interrupted  by  the  phasinq  sector  for  an  angle  of  2  arcsin  r(=2ft) 
for  each  rotation  of  the  reticle.  For  this  particular  fiqure  only,  the 
imaqe  has  been  rotated  around  the  centre  of  the  stationary  reticle. 


FIGURE  3  -  Reticle  geometry 
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Suppose  that  the  image  is  at  P(r,'i  ,  shown  in  Fin.  4(a),  which 
has  been  drawn  for  the  instant  that  the  retielp  and  spatial  references 
are  coincident.  The  photocell  sional  will  then  be  as  in  Fin.  4(b) 
Reference  signals  (Fiq.  4(c))  are  generated  by  the  reticle  motor,  and 
are  phased  so  that  the  information  contained  in  the  photocell  signal  can 
be  correctly  recovered.  The  photocell  signal  is  rectified  (Fig.  4(d)) 
and,  in  conjunction  with  the  reference  signals,  its  SINE  and  COSINE 
components  are  determined. 


1 

“  2 IT-  $-4; 

2r  1 

COSINE  =  — 

feus  a  da  + 

1  cos  a  dot 

TT 

J 

L  o 

Tv+Q-P  J 

7 

7T 

sm  6  cos  ip  , 

2 

=  -  —  r  cos  ip  , 

IT 

where 

a  =  to  t  . 
s 

Similarly 

SINE  =  —  r  sin  .  [2] 

IT 

The  spatial  axes  are  defined  such  that  when  ^  =  0,  the  imane  lies  on  the 
positive  X  axis,  and  when  =  w/2,  it  lies  on  the  positive  V  axis.  It 
is  required  that,  when  an  error  is  in  the  positive  direction,  the  error 
signal  be  also  positive.  For  this  to  occur,  the  error  siqnals  for  the 
X  and  Y  directions  (E^  and  E  )  must  be  equal  to  (-COSINE)  and  SINE 
respectively,  which  requires  the  use  of  (-cosine)  and  sine  reference 
signals.  Henceforth,  these  will  be  the  reference  sinnals  employed. 
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TARGET 

IMAGE 


FIGURE  A  -  Waveforms  required  for  example 
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As  this  simplified  example  demonstrates  a  linear  control  it  would 
be  anticipated  that  the  more  sophisticated  model,  which  is  emnloyed 
throughout  the  remaining  analyses,  will  show  a  similar  tendency. 


2.3  Signal  Processing 


The  photocell  signal  is  first  filtered.  In  the  model,  only  the 

carrier  (at  the  spoke  frequency)  and  its  immediate  side-bands  are  nassed 

by  the  filter.  If  the  spin  frequency  is  fg,  that  of  the  carrier  will  be 

12f  (assuming  a  12— spoke  reticle),  and  the  filter  will  pass  only 
s 

components  at  llfs>  12f  and  1  -5f ^ .  In  oractice,  the  sianals  from  the 
various  target  images  are  summed  by  the  photocell,  and  then  filtered. 
As  the  signal  processing  to  this  point  is  linear,  the  model  achieves  the 
same  result  by  filtering  the  siqnals  from  the  individual  images,  and 
then  summing.  Thus  the  filtered  siqnal  from  the  i'th  inane  ran  be 
represented  by 


13 

S(F . )  =  l  (a.,  sin  ja  +  b. .  cos  ja)  • 
1  j=ll  ^  XJ 


[3] 


The  total  filtered  siqnal  is  then 

S (F)  =  I  S(F  )  .  [4] 

i 

Before  this  siqnal  is  detected,  some  form  of  automatic  oain 
control  must  be  applied,  otherwise  the  amplitude  of  the  error  sinnals 
will  be  a  function  of  the  imaqe  intensities,  and  linear  control  would 
not  be  feasible.  The  type  of  AGC  employed  herein,  is  that  nf  dividinq 
S(F)  by  its  root-mean-sauare  (RMS)  value,  to  give 


S(AGC)  =  S(F)/RMS, 


[5] 
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where 

13  2  2  * 

RMS  =  [  £  (a  +  b‘)]  ,  r6l 

j=ll  3  3 

and  a.  =  I  aij  ;  b.  -  Z  b..  . 

This  sinnal  is  then  square-law  detected,  and  the  coefficients  of  sin  i 
and  cos  a  extracted.  Let 


13 

S(AGC)  =  £  (c .  sin  ja  +  d.  cos  ja)  .  [7] 

j=ll  3  3 

Then 

S2(AGC)  =  Ex  cos  a  +  Ey  sin  a  +  higher  frequency  components,  [8] 


where 

EX  =  "(c12tcll  +  C13)  +  d12(dIl  +  d13^  ' 
EY  =  C12(dll  "  d13)  ”  d12(Cll  ‘  C13}’ 


The  higher  frequency  conoonents  are  filtered  out,  and  Cy  and  Fy  ire  used 
to  drive  the  reticle  assembly,  in  the  X  and  Y  directions  resnect ivelv , 
via  the  tor cue  motors  which  have  a  gain  of  C. 


A  flowchart  of  the  model  is  given  in  Fig.  5. 

2.4  Error  Sional__fo£  Small  Error 

If  the  target  image  is  displaced  only  a  small  distance  from  the 
reticle  centre,  then  sine  and  cosine  approximations  can  be  emnloyed  to 
simplify  the  calculations .  This  will  permit  a  further  examnle  of  the 
signal  orocessinq  to  be  qiven,  without  as  yet  resortino  to  a  computer. 
In  Fig.  6(a)  is  shown  the  photocell  signal  for  a  target  imane  on  the  X 


n 
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FIGURE  b  -  Flowchart  of  seeker  model 
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E  sin(n  +  p)  a  +  cos (n  +  p)a. 


The  waveform  of  Fig.  6(a)  will  be  analysed  by  first  separatinq  it  into 
the  two  parts  shown  in  Fiq.  6(b).  It  follows  that  the  S(F)  for  the 
photocell  signal  is 


S(F) 


1  -  cos(p  -  l)c 

(P  -  Du 
0 


-2  1  -  cos  pB 

—  + - 

it  pir 

0 


1  -  cos(p+l)B 

(P  +  Dir 
0 


[11] 


Assume  nr<<l,  so  that  B  =  (arccin  r)  *  r,  then  eo.  11  can  be  written 


S(F) 


(p  -  1)  r2 
2ir 
0 


-2  (p  +  1)  r2 

n  2n 

0  0 


[12] 


Let  the  target  image  be  displaced  from  C(x,0),  or  P(r,0),  to  C(x,y),  or 
P  ( r  ,41 ) ,  where 

r  =  (x2  +  y2)*  , 
ip  =  Arctan  y/x  . 


In  the  model  the  phase  angle  is  introduced  after  the  filtering,  so  that 


2 

(p-1)  r  cos(p-l)* 
2  ir 


S(F)  = 


(P-D  r  sin(p-l)sp 
2  ir 


-2^  cos  pip 

it 


-2_  sin  pif/ 
ir 


2 

(p+1)  r  cos  (p+1) 'J' 

2  n 


2 

(p+1)  r  sin(p+l)i// 
2  v 


In  the  case  of  multiple  tarqets,  it  is  at  this  ooint  that  the  S(F.) 
would  be  summed  to  form  S(F).  The  RMS  value,  because  r  is  small,  is 
very  close  to  2/it .  Hence 
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S(AGC)  = 


2  2 

(p-1)  r  cos(p-l)iJ/  -cos  pif i  (p+1)  r  cos(p+l)i|* 

4  4 

2  2 

(p-1)  r  sin(p-l)i|/  -sin  pi|i  (p+1)  r  sin(p+l)tj/ 

4  4 


[14] 


Equation  14  is  substituted  in  ea.  8  to  give 


EX  = 


rz  {cos  pip [  (p-l)cos(p-l)ij>  +  (p+l)cos(p+l)\j/] 

4 

+  sin  pip  [  (p-l)sin(p-l)\(  +  (p+l)sin(p+l)ij>] } 


[15] 


=  6  r  cosip 


E  _  {cos  p^  [  (p-l)sin(p-l)^  -  (p+l)sin(p+l)\p] 

Y  4 

-  sin  pip  [  (p-l)cos(p-l)*  -  (p+licosfp+l)^]} 
2 

=  6  r  sin<|/ 


[16] 


It  is  seen  that  for  small  disDlacements  the  system  has  a  sciuare-law 
control . 

2 . 5  Equilibrium  Position  for  Two  Targets  (Small  Displacement) 


A  further  example  is  now  qiven  where  two  tarqets  are  oresent,  but 
with  only  a  small  separation  between  them.  Let  the  first  tarnet  have  an 
intensity  of  unity  and  a  position  at  eouilibrium  of  C(x  ,n),  and  the 
other  target  an  intensity  of  N  and  a  nosition  at  eouilibrium  of 
C(-X2»0).  Then  if  pxj(x2)<<;l»  hv  annlyinq  eo.  12  to  both  tarnpts,  the 
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filtered  photocell  signal  is 


S(F) 


(P-1) (x2-Nx2)  -2(N+1) 

2tt  it 

0  0 


(Eii)  (x2-Nx2) 
2m 


0 


For  eauilibrium  must  be  zero 


(refer  to  eo.  8),  that  is 


2(N+l)p(x, 2  -  NX  \  =  0, 

2 

IT 

the  solution  of  which  is 


[17] 


[18] 


[19] 


This  is  to  be  expected  from  the  sauare-law  relationship  between  the 
error  signal  and  the  image  error  (r)  derived  in  Section  2. A.  This 
position  can  be  seen  to  be  stable,  as  eo.  18  shows  that  an  increase  of 
Xp  with  a  conseguent  decrease  in  x^>  will  cause  E^  to  become  positive, 
whereas  if  x^  becomes  negative  so  does  E^.  A  positive  error  signal  will 
cause  the  reticle  to  move  to  the  right  and  hence  decrease  the  error  (x^) 
between  the  reticle  centre  and  a  taroet  image  lying  to  the  right  of 
centre . 


t 

i 

i 
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2 . 6  txtended  Target 

First  consider  the  case  where  two  tarqets,  at  the  same  distance 
from  the  centre,  are  separated  by  the  width  of  one  sector,  as  shown  in 
Fig.  7(a),  which  has  been  illustrated  to  show  the  instant  at  which  the 
reticle  and  spatial  references  are  coincident.  The  waveforms  resulting 
from  the  two  individual  images  are  shown  in  Fig.  7(b),  and  the  total 

signal  in  Fig.  7(c),  where  and  $g  are  the  arqunents  of  the  nositions 

of  images  A  and  B. 

An  analysis  of  the  waveform  of  Fia.  7(c),  which  is  given  in 
Appendix  A,  shows  that 

2it  2ir 

sin  pa  •  f(A+B)  da  =  ^  cos  pa  •  f(A+B)  da  =  0  .  [20] 

and  thus  a^  =  bu  =  0;  it  can  be  seen,  from  eo.  8,  that  because  of  this 

the  error  signal  will  also  be  zero.  It  follows  that  if  the  target  image 
was  in  the  form  of  an  arc,  covering  more  than  one  sector,  the  carrit  ' 
signals  from  all  points  separated  by  exactly  one  sector's  width  would 
nullify  one  another.  If  the  arc  was  to  cover  exactly  an  even  number  of 
sectors,  the  carrier  would  be  lost,  and  the  error  signal  would  be  zero. 

The  signals  from  an  arc  image  will  now  be  examined.  An  arc  is 
chosen  as  an  image  shape,  in  preference  say  to  a  straight  line,  onlv 
because  of  its  mathematical  simplicity  in  this  introductory  analysis. 
It  suffers  from  the  disadvantage  that  nothinq  can  be  concluded  about  the 
stability  of  equilibrium  positions,  as  once  the  image  moves  it  is  no 
longer  an  arc  with  the  same  centre,  and  hence  the  eouations  emploved  in 
the  analysis  become  invalid. 
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FIGURE  7  -  Waveforms  from  two  images 
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In  Fin.  8(a),  the  taroet  imaqe  is  shown  by  a  heavy  are  from 
P(1,'A)  to  P(1,A).  To  derive  the  error  sionals  resultinn  from  this 
imaqe,  the  filtered  siqnal  for  an  imaqe  of  unit  intensity  situated  at 
P(1,0)  will  first  be  calculated.  Its  nhotocell  siqnal  is  shown  in 
Fig.  8(b),  and  hence  (see  also  Appendix  B) 


S(F) 


u 


11 

0 


u 


12 

0 


u 


13 

0 


1 2i ; 


where 

un  =  7.596/llir, 

U12  = 

u13  =  7 . 596/ 13m . 

If  the  intensity  of  the  arc  is  uniform,  and  its  total  radiance  is  unity, 
then  for  a  segment  at  P ( 1 , X )  of  Ipnqth  6X 


SCF6X> 


<SX 

Uj^COS 

11X 

u^cos  12X 

U13' 

2A 

, 

u^sin 

11X 

u^sin  12X 

U13 

[22] 


In  the  limit 


S(F) 


ii 

/  S(F^)  as  6A-rdX->-0, 


1_ 

A 


U11  sin  11A  J^12  sin  12A  U13  sin  13A 

11  12  13 


[23] 


0 


0 


0 
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Ihe  carrier  amplitude  of  the  arc  is 

=  (sin  12A)/12itA,  [24] 

and  the  sum  of  the  side -band  amplitudes  is 

=  7 .596  ( sin  11A 

S  -  rrA  \  n2 

These  arp  plotted  in  Fin.  9.  As  the  peaks  of  both  the  carrier  and  the 
sum  of  the  side-bands  decrease  raoidlv  with  arc  length,  their 
contribution  to  the  total  error  signal,  in  the  presence  of  a  point 
source  of  similar  intensity,  will  be  very  small. 

Further,  an  arc  of  total  length  sav  n  would  be  in  enuilihrium,  as 
the  carrier  amplitude  is  then  zero.  A  point,  image  at  the  centre  of  the 
reticle  is  also  in  eguilibrium,  and  eo.  12  shows  that  this  is  because 
the  side-band  amplitudes  are  zero.  However,  if  both  the  arc  and  point 
images  are  present  at  the  same  time,  the  seeker  is  no  longer  in 
equilibrium,  as  the  side-bands  from  the  arc  and  the  carrier  from  the 
point  will  generate  an  error  signal  and  cause  a  displacement  of  the 


sin  13a\ 
132  / 


seeker . 
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3.0  RESULTS 


3.1  General 


The  preceding  analyses  have  served  either  to  illustrate  the 
functionino  of  the  seeker  or  to  provide  a  means  of  validating  the 
computer  orogramme,  and  as  such  have  investigated  special  cases  or  used 
approximations.  The  rest  of  the  study,  which  will  be  qeneral,  will 
reauire  the  use  of  a  diaital  computer.  The  simulation  model  was 
programmed  in  APL,  and  run  on  the  Xerox  560  digital  computer.  The 
programmes  are  listed  in  Appendix  C. 

3.2  Error  Curves 


The  error  curve  (Fig.  10(a))  gives  the  error  signal  as  a  function 
of  the  distance  of  a  point  imaae  from  the  centre  of  the  reticle.  It  is 
to  be  noted  that  the  curve  gives  an  overall  aopearance  of  linearitv,  as 
was  suggested  by  the  reasonina  in  Section  2.2.  The  number  of  neaks  on 
the  curve  corresponds  to  the  number  of  changes  of  snoke  interceptions 
with  the  phasing  sector  as  r  varies  from  0  to  1. 

Fiqure  10(b),  which  is  plotted  from  the  data  of  Table  I,  shows 
the  portion  of  the  error  curve  for  small  values  of  r.  In  Section  2. A 
the  error  signal  was  derived  as  being  eaual  to  6r^  for  values  of  r  less 
than  0.08;  this  verifies  that  the  Dart  of  the  error  curve  for  small 
values  of  r  is  correct. 


i 
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TABLE  I 

Error  signal  for  small  r 


m 

X 

6r2 

01 

0.000599 

0 . 0006 

02 

0.00239 

0 .0024 

03 

0.00536 

0.0054 

04 

0.00946 

0.0096 

OS 

0.0147 

0.015 

06 

0.0209 

0.0216 

07 

0.0281 

0 .0294 

03 

0.0362 

0.0384 

09 

0.045 

0.0486 

1 

0.0545 

0 . 06 

11 

0.0646 

0.0726 

12 

0.075 

0.0864 

13 

0.0857 

0.101 

14 

0.0965 

0.118 

IS 

0.107 

0.135 
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5.3  Eciuil ibrium  Position  for  Two  Point  Targets 


In  this  test,  two  point  target  imaqes  are  nlaced  at  a  oiven 
distance  apart  and,  after  equilibrium  has  been  reached,  their  distances 
from  the  reticle  centre  are  noted.  The  two  imaqes  are  situated  at 
CCx^jO)  and  C(x2*0),  with  intensities  of  1  and  N.  The  imaqes  are  a 
fixed  distance  apart  (d),  where  d  =  x^  -  x2- 


From  Fig.  11  and  Table  II  it  can  be  seen  for  rf  =  0.001  that 


or 


log  N  |x2/x1|  =  i  log  N, 

jx1/x2(  =  N*. 


[26] 


This  result  is  in  agreement  with  eo.  19. 


For  larger  d  the  value  of  N|  x7/xjJ  is  rlose  to  unitv,  that  is 
|xj/x7|*N,  if  neither  image  is  close  to  the  centre  of  the  reticle,  so 
that  a  weighted  mean  can  be  applied  in  this  case  to  calculate  the 
equilibrium  position.  Once  either  of  the  imaoes  is  close  to  the  reticle 
centre,  the  signals  from  this  imaqe  will  no  lonner  be  linearlv 
proportional  to  its  distance  from  the  centre;  nevertheless,  the  slones 
of  the  curves  for  large  N  remains  one  half. 


UNCLASSIFIED 

26 


LEGEND  .. 


10  r 


5 


3 

2 


—  d  =  0.001 

—  d  =  0.4 
._  d  =  1.0 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


X 


X 


X 


X 


X 


/ 

/ 

s 

/ 


/ 


x 


1  . .-.**£.1 . . 

I  10  100 

N 


FIGURE  11  -  Two  point  targets 
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TABLE  II 


Iwo  ooint  imanes 


N 

_ 

X 

2 

jx  /X  ] 

1  j/  2I 

N|yx2i 

d  =  0.001 

1.00 

.00050 

-  .00050 

1 . 00 

1 . 00 

1.58 

.00056 

-  .00044 

1  .26 

1  .  26 

2.  SI 

.00061 

-  .  00039 

1 .58 

1  .59 

3 .23 

.  00067 

-.00033 

1  .  99 

2 . 00 

6.31 

.00072 

-  .00028 

2.51 

2.51 

10.03 

. 00076 

-  .00024 

3.16 

3.16 

15.8$ 

.  00080 

-  .00020 

3.98 

3.98 

25.12 

.03333 

-.00017 

5.01 

5.01 

39.81 

.00036 

-  .00014 

6.31 

6.31 

63.10 

.00039 

-  .0001  1 

7 . 94 

7 . 94 

100.00 

.  00091 

-  .00009 

10.00 

10.00 

d  =  0.4 

.1.00 

.  200 

-  .200 

1 . 00 

1 . 00 

1  .  26 

.  225 

-.175 

1 .28 

.93 

•  1.58 

.  249 

-  .151 

1  .  64 

.  97 

2.00 

.  269 

-.131 

2.07 

.  97 

2.51 

.  235 

-  .  1 1  5 

2.48 

1.01 

3.  16 

.  293 

-  .  102 

2.92 

1.08 

■  3.93 

.  309 

-  .  031 

1  .  00 

3 . 93 

6.01 

.318 

-  .082 

3.83 

1  .  29 

6.3  1 

.  326 

-  .  074 

4 . 44 

1.42 

10.00 

.34  0 

-  .060 

5.72 

1.75 

15.85 

.352 

-  .04  3 

7.3  3 

2.16 

25.12 

.  36  1 

-  .039 

9.31 

2.70 

39.81 

.36  9 

-  .031 

11.78 

3.38 

6  3.10 

.  375 

-  .025 

15.00 

4.21 

100.00 

.380 

-  .  020 

18.80 

5.32 

d  =  1.0 

1.00 

.500 

-  .500 

1.00 

1.00 

1  .  26 

.  576 

-  .424 

1.36 

.  93 

1  .  5  3 

.613 

-.337 

1  .  59 

1  .  00 

2.0  0 

.  649 

-  .351 

1.85 

1.03 

2.51 

.  693 

-  .  307 

2  .  26 

1.11 

3.  16 

.  781 

-.21? 

3 . 56 

.  89 

3.93 

.817 

-  .  183 

4.45 

.  89 

5.01 

.843 

-  .  157 

5.37 

.  93 

6.31 

.  865 

-.135 

5.41 

.  98 

10.00 

.  897 

-  .  103 

8 . 67 

1.15 

15.85 

.  919 

-  .081 

11.29 

1.40 

25.12 

.915 

-  .065 

14.48 

1.73 

39.81 

.  949 

-  .051 

18.49 

2.15 

63.10 

.959 

-.041 

23.57 

2.68 

100.00 

.96  3 

-.032 

30.06 

3.33 
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3 . 4  Representation  of  Line  Targets 


For  various  numbers  of  points  formino  the  line-tarnet  imaqe,  the 
values  of  a^,  a|2  anc*  a]3  are  calculated  and  oresented  in  Table  III  and 
Fig  .  12,  where  the  line  is  0.8  long  and  positioned  symmetrically  about 
the  X  axis  at  x  =  0.5.  It  was  concluded  from  the  figure  that  little 
error  should  occur  by  employing  21  points. 


To  verify  that  this  number  nf  points  is  sufficient  to  renresent  a 
line  image  anywhere  on  the  X  axis,  the  amplitudes  of  the  carrier  and  the 
sum  of  the  side-band  amplitudes  were  plotted  for  the  same  line  imaoe, 
and  its  position  along  the  X  axis  was  varied  from  -1.0  to  1.0.  The 
graph  is  given  in  Fig.  13  for  a  line  image  represented  by  21  ooints  and 
in  Fig.  14  for  41  points.  The  only  discernable  difference  between  the 
two  curves  is  in  the  peaks  of  the  carrier  immediately  on  either  side  of 
the  ordinate. 


TABLE  III 

Number  of  points  to  form  a  line  image 
Image  length  =0.8  (Symmetrical  about  X  axis), 


Position  of  image:  x  =  0.5  . 


NUMBER  OF 
POINTS 

ail 

S12 

a13 

1 

.1099 

-  .  5305 

.  0929 

3 

.0816 

-  .4953 

.1065 

7 

.0232 

-.0882 

.0113 

11 

.0219 

-.0821 

.0100 

21 

.0213 

-.0/92 

.0094 

31 

.0212 

-.0786 

.  0093 

51 

.0211 

-.0783 

.0093 

101 

.0211 

-.U7H2 

.0092 
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FIGURE  13  -  Line  signals  for  21  points 


1  -1 


FIGURE  14  -  Line  signals  for  41  points 
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The  difference  between  the  signals  from  the  representation  of  a 
straight-line  image  by  21  and  41  points,  as  shown  in  Figs.  13  and  14,  is 
small,  and  hence  the  following  computations  have  been  done  with  21 
points;  the  effect  of  the  difference,  if  any,  will  be  discussed  at  the 
end  of  this  section. 

In  Figs.  15(a)-17(a)  are  shown  the  values  t,f  a^  and  ag(=  an  + 

)  for  a  point  image  of  intensity  N  and  a  line  imane  of  unit  intensity 
and  of  length  0.8  centred  about  the  X  axis.  It  is  now  required  to 
derive  the  points  of  equilibrium  from  these  curves  for  a  given 
separation  of  the  two  imaqes.  The  separation  will  be  taken  to  be  0.5, 
with  the  point  image  lying  to  the  left  of  the  line  imane.  For 
equilibrium  to  occur  the  error  signal  must  be  zero,  that  is 

[al2(xp)  +  ai2(xL)]  tas(xp)  +  as^xL^  =  0  ’  '•27^ 


where 


Xp  is  the  position  of  the  point  image, 

x^  is  the  position  of  the  line  inane, 

a-(x.  )  is  the  value  of  a  •  at  the  point  x, 
J  k  J  K 


on  the  k  curve, 


and 

xL  -  xp  =  0.5  . 


ti 
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FIGURE  15  -  Line  and  point  images 
separation  =  0.5,  N  = 
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- -13 


-  .751 


LEGEND. 
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FIGURE  17  -  Line  and  point  images 
separation  =  0.5,  N  = 
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Thus  the  equilibrium  position  can  be  given  by  the  solutions  of  both 


a!2(xL  '  0l5)  +  a12(xL3  =  0  * 
ag(xL  -  O.S)  +  ag(xL)  =  0  . 


[28] 


These  equations  could  he  solved  bv  plotting  the  error  siqnal,  and 
noting  the  zero  crossings.  A  positive  slope  at  these  positions  would 
show  stability,  as  a  movement  of  the  reticle  centre  to  the  left  would  be 
accompanied  by  a  positive  error  siqnal,  which  would  tend  to  return  the 
reticle  to  its  equilibrium  position.  Althouqh  the  error  siqnal  can  he 
very  easily  calculated,  the  value  of  RMS  is  involved,  for  which  a^,  ap 
and  a^  are  required.  A  more  direct  approach  is  to  use  the  data  fron 
the  curves  of  Figs.  15(a)-17(a),  where  only  the  a^'s  and  ss's  are 
given.  Let 

Ex  =  E^/RMS2  ,  [29] 

2  1 

Then  because  RMS  cannot  be  negative,  E^  must  have  the  samp  zero 

crossings  as  Ex,  and  the  signs  of  the  slnnes  at  these  positions  must  be 

the  same  for  rv  and  f'  Thus,  the  eouilibrium  positions  are  determined 
A  X 

I 

by  plotting  E^,  where 

EX  =  [al2(XL  '  °'5)  +  al2(xL)]  USCXL  -  °’53  +  aS(xL3]  ‘  [30] 


In  Figs.  15(b),  16(b)  and  17(b)  this  curve  is  shown  for  N  =  l.D,  0.5  and 
0.1  respectively. 

It  can  be  seen  for  N  =  1.0  that  there  is  one  eouilibrium 
position,  that  it  is  stable,  and  occurs  for  =  0.A5,  Xp  =  -0.07.  Fnr 
N  =  0.3,  there  are  three  equilibrium  positions,  two  nf  which  are  stahle. 
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For  N  =  0.1,  there  are  15  equilibrium  positions  of  which  8  are  stable. 
The  figures  show  that  if  a  41-noint  representation  had  been  employed 
instead  of  a  21-point  one,  the  larger  carrier  peak  to  the  right  of  the 
ordinate  might  have  caused  two  more  equilibrium  positions  to  have 
occurred  at  x^s  between  0  and  0.1,  for  N  =  0.3. 
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A .  0  DISCUSSION 

In  practice,  a  small  nortian  of  the  centre  of  the  reticle  may  be 
blacked  out  to  prevent  oscillation  of  the  seeker  when  the  error  is 
small.  This  would  change  the  shape  of  the  two  point  target  eauilibrium 
curves  in  the  region  where  one  of  the  imaaes  enters  this  zero 
transmission  zone. 

Multiple  statically  stable  equilibrium  positions  were  found  under 
certain  conditions  when  both  a  point  and  a  line  imane  were  present  on 
the  reticle.  Under  dynamic  conditions  some  of  these  positions  may  not 
be  stable. 


The  filters  employed  extracted  the  carrier  and  its  immediate 
side-bands  with  no  attenuation  or  phase  shift.  The  simulation  model 
could  be  readily  adapted  to  filters  requiring  attenuation  and  phase 
shift  at  particular  harmonics. 

The  only  extended  target  examined  was  a  line.  However,  areal 
tarqets,  such  as  illustrated  in  Fiq.  18,  could  be  studied,  but  because 
of  the  large  number  of  points  required  to  define  an  areal  target  and  the 
conseauent  large  number  of  computer  operations  entailed,  the  programming 
would  probably  have  to  be  performed  in  a  language  other  than  APL. 

The  reticle  employed  is  the  'sunburst'  reticle  with  straight 
spokes.  Other  reticles  could  be  employed  with  very  small  changes  to  the 
simulation.  One  such  reticle  is  illustrated  in  Fiq.  19,  in  which  the 


spiral  spokes  improve  the  spatial  filtering. 
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5.0  CONCLUSION 


A  simulation  model  of  a  generic  seeker  employing  an  AM  reticle 
was  developed,  from  which  a  digital  programme  was  prepared  to  calculate 
the  seeker  response  to  various  stimuli.  It  was  demonstrated  that  the 
eguilibrium  position  of  two  point  images  can  be  described  by  a  weighted 
mean  only  if  the  separation  and  intensity  ratio  is  such  that  neither 
image,  at  equilibrium,  is  nearer  the  reticle  centre  than  0.1  of  the 
reticle  radius.  This  occurs  because  the  system  is  linear  except  for 
images  near  the  centre  of  the  reticle. 

It  was  shewn  that,  because  of  phase  cancellations,  an  extended 
image  contributes  less  to  the  error  signal  than  a  point  of  the  same 
total  intensity;  hence  a  weiahted  mean  cannot  be  used  for  determininn 
the  eguilibrium  position  of  an  extended  and  a  ooint  imane.  The 
programme  disclosed  that  with  both  a  line  and  a  point  imaoe  nresent  on 
the  reticle,  for  certain  ranges  of  intensity  ratio,  there  could  he  many 
static  equilibrium  positions. 

The  model  could  be  modified  to  accommodate  the  given 
characteristics  of  a  particular  AM  reticle  seeker,  such  as  the  neometry 
of  the  reticle,  the  amplitude  and  phase  responses  of  the  carrier  filter, 
and  the  design  of  the  automatic  gain  control.  Areal  targets  of 
specified  shapes  and  gradations  in  intensity  could  also  be  simulated. 
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APPENDIX  A 


From  Fig.  7 


Derivation  of  Equation  20 


f(B,a)  =  f(A,  a+n/12'), 


[Al] 


where 


a  =  to  t  . 
s 

It  follows  that 


2ir  2  w 

.X*  f(B,a)sin  12a  da  =  J"  f(A, 


a+n/12)sin  12a  da. 


2ir 


-  J  f  {A,a)sin(12a  -  n)  da, 
*0 


[A2] 


2ir 


--J  f(A,a)sin  12a  da. 
0 


Theref  ore 


2it  2it 

f(A+B)sin  12a  da  =  /  [f(A)  +  f(B)]sin  12a 

r\ 


da 


[A3] 


Similarly 


2ir 

y"  f  (A+B)  cos 
0 


12a  da  *  0. 


[A4] 
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APPENDIX  B 


Derivation  of  the  Constant  in  Equation  21 

From  Fiq.  8 


1_ 

np 


2 

2 

irp 


(photocell  signal)sin  pa  da. 


I  C-l)i[cos(Zi^P7r  -  cos*-j^ptt] 
i=0  1 *  11 


(1-cos 

(1-cos 


l  cos(7+2j)^y  , 
j=0 

H|)cos  (1+2cos|2IL). 


Therefore 


a12  =  7.596/llir, 

a12  =  '  1/ir* 
aJ3  =  7.596/13*. 


[Bl] 


[B2] 
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APPENDIX  C 


Listing  of  Computer  Programmes 


Function  QL2 


VOUT-QL2  ;I;  C;  /I;  P;  5;  T;  RMS 

[1]  I- 1 

[2]  P«-C’-0 

[33  wc:r-< <mr:t*A'x[r])+mn*m.n )*.5 

[4]  -*-(2’>i  )//-'r 

[s;i  ozz[r]*cti  xxcn.rym 

[6]  A<-1  2 °.oGC*FI 

[7]  P-P+(G[l;]*/>[2;])-C[2i]*.A[l;] 

[a]  ;*•$+(£[  1{]M[1;])+C[2;]M[2;] 

CO]  V2’:-(7  =  pAA)/PA 

CIO]  /«-/+•  1 

[11]  -*.V£ 

[12]  P.V:P-*f5-(+/Px5^P.i)*.  5 

[13]  P*- Pi  RMS 

[14]  U+SiRMS 

[15]  P-+/"l  +  Pxl4>P 

[16]  S+S*+/~  1U,>1$2 

[17]  2'*-+/”l  l4*P 

[18]  7«-2'-*/-l*P*l<K 

[19]  OUT+-(. -S)  .7 
V 


Use:  Calculates  E^  and  Ey  for  any  number  of  point  images  defined  by  the 
vector  XX,  YY  and  ZZ,  where  XX(I),  YY(I)  and  ZZ(I)  correspond  to 
the  position  and  intensity  of  the  i'th  image. 


Calls:  QL,  QL1 


UNCLASSIFIED 

A3 


Function  QL 1 


RG<-QL\  RRiYiPiQiCOiEiTERViN 
[  1  ]  Y+QL  RR 

[2]  F*-y[l;] 

C3]  a«-y[2;] 

[4]  C0+1OP[1]xCG-(  .  5*^24)  +  (  -1+L  .  5  *  TERM )  +  \  TERM*-  3 

c  s ]  y-(o.c[i]=o)** 

[6]  P«-r[l;] 

[7]  £>-$[l]$((L.5xpP) ,2)pP 

[8]  G«-+/[2]-/[2]2  1*  .OE • .  *GG 

C  93  G[2;]«-G[2;]+C0 

CIO]  <?[  1  j  ]*--<7[  1 ;  ] 

[11]  G*-Cto(2,TERH)  pGG 
V 


Use:  Determines  a^,  (i  =  11,  12,  13)  for  a  given  point  image 


Calls:  QL1 


Called  by:  QL2 


UNCLASSIFIED 

44 


Function  QL 


V'l*-QL  R a ;  h ;  A  ;  7. ;  P ;  C ;  £  ;v ;  H 

[  1  ]  /?•*•  (  +  /  R  a  *  R  A  )  *  .  5 

[  2  ]  FI-H  ft  [  1  ]  /  T  A  i'J  HR  [  2  ] 

[3]  A«-  (  Of  .  5  *  ft )  *  0  ,  i  ~  1  +  //*■  2  4 

C  4  ]  X*-  (  1 ,  pA  )  pA 

L  5  3  A*- A  ,  [  1  ]  (  ,  pA )  pO  1 

C  6  ]  A'C  1 ;  ]-‘-A[  l;]+(o2)*  P-*-A[  1 ;  ]  <0 
[  7  ]  A+-  (  +  /  P)  i>A 

[8]  A[l;]~A[l;]-(o2)*P*-A[l;]>o2 

C9]  Y«-(-  +  /F)*  A 

Cio]  y*-(o,  ( ( i+f>y)  pi )  ,o)\y 

[n]  Y[l;l>0 

[12]  Y[l;l+pY]-o2 

[13]  Y L  2  ;  1  ]«-l  -  Y [  2  ;  2  ] 

[14]  Y[2;l*pY]-l-Y[2;“l+lfpY] 

[15]  2-<c.5)-_2 OP 

[16]  F-*-Y[  1 ;  ] 

[17]  £-P>2: 

[18]  ife«-Cil 

[19]  -(fc^=l)//10 

[20]  CCC£~1  ]«-l 

[21]  y^t/y 

[22]  y[l;l>Z 

[23]  AD:2-(02)-Z 

[24]  F*- y [  1 ;  ] 

[25]  w-P<Z 

[26]  (i'<*C\Q 

[27]  -($i)>pS)/0 

[28]  £[fi5>l 

[29]  Y-C/Y 

[30]  Y[l;“l*pY]«-2 
V 


Use:  Determines  zero-crossing  positions  of  the  photocell  waveform  for 

a  given  uoint  image . 


Called  by:  QL1 


UNCLASSIFIED 

45 


Function  MP 


VI«-.VP  NiAiJiMOliPOl 
[1]  M0\+lItiE\.',2  3] 

C  2]  P01«-PIiY[;2  3] 

C  3  ]  I*  1 

[  4  3  A*\ 0 

[5]  4S:<7>*I>50 

[6]  A*- A  ,  (A/01  [I;  1  3  +  P0l[<7;  1  3  *«)  *  A/01  [  I ;  2  3+ P01  [  </;  2  3* // 

[73  -*(<7=200  )/d0 

[83  1*1*1 

[93  -*>15 

[103  ADiI<—A 
V 


Use:  Calculates  the  data  for  Fiqs  15(b),  16(b)  and  17(b)  from  the 

matrices  LINE  and  PIN,  where  LINE  contains  the  values  of  a^x^ 

and  a<j(x^)  for  -lsx^sl,  and  PIN  contains  the  values  of  a^Cx^) 

and  a  (x  )  for  -lSx  si. 

5  P  p 


UNCLASSIFIED 

46 


Function  ATAN 


VZ-X  ATAN  Y 

[1]  ■*(  IE~6<\X)/U1 

[ 2 ]  X-0 

[3]  Ul:Z+-~3oY*X+lE  7 

[4]  -u<o)/yp 

C5]  +AB 

[6]  yP:Z*-Z+01 

[7]  AB:Z->-{02)\Z 

[8]  ■♦(  Zso  1 )  /0 

[9]  Z-Z-02 
V 


Use:  Calculates  the  arctangent  and  assigns  to  it  the  appropriate  value 

dependent  on  its  quadrant. 


